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Abstract 

O 

CN . 

The role of anharmonic effects on the excitation of the double giant dipole 

CN 

resonance is investigated in a simple macroscopic model. Perturbation theory 
is used to find energies and wavefunctions of the anharmonic oscillator. The 
cross sections for the electromagnetic excitation of the one- and two-phonon 

a 

giant dipole resonances in energetic heavy ion collisions are then evaluated 
through a semiclassical coupled-channel calculation. It is argued that the 
variations of the strength of the anharmonic potential should be combined 
with appropriate changes in the oscillator frequency, in order to keep the 
giant dipole resonance energy consistent with the experimental value. When 
this is taken into account, the effects of anharmonicities on the double giant 
dipole resonance excitation probabilities are small and cannot account for the 
well known discrepancy between theory and experiment. 
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The double giant dipole resonance (DGDR) has attracted considerable interest in the last 
decade. Several experiments to measure the DGDR cross section using relativistic heavy 
ion beams have been performed [0-^. Comparison with the predictions of the harmonic 
oscillator model has clearly demonstrated a systematic discrepancy. The experimental values 
for the DGDR cross sections exceed the theoretical predictions by a considerable amount. 
One of the attempts to explain these differences was made by Bortignon and Dasso [0, 
using a macroscopic anharmonic oscillator model. These authors found that with a small 
anharmonic perturbation of the r 4 -type one can reproduce both the experimentally observed 
DGDR excitation energy (which only marginally differs from that obtained in the harmonic 
approximation) and the DGDR cross section for the 208 Pb + 208 Pb collision at 640 A-MeV. 
They reached a similar conclusion for the 136 Xe + 208 Pb collision at 700 A-MeV, where a 
much greater discrepancy from the harmonic model appears 0. The purpose of this paper 
is to point out that this model does not lead to the enhancement found in Ref. if 
proper renormalization of the oscillator frequency is performed in order to guarantee that 
the theoretical giant dipole resonance (GDR) excitation energy is kept at the experimental 
value. 

The model of Refs. |7],?] is based on the following Hamiltonian 

H = H + F{x,y,z;t), (1) 
where H is the anharmonic oscillator describing the intrinsic motion of the projectile, 

Ho = ± (pI +pI+pI) + ^ {x 2 + y 2 + z 2 ) + f {x 2 + y 2 + zY, (2) 

where D is the mass parameter, C is the oscillator strength and B is the strength of the 
anharmonicity. Here, we take the mass parameter to be the reduced mass for the motion of 
the protons against the neutrons, 

NZ 
= ~A~ m0) 

where is the average nucleon mass. The beam is assumed to be parallel to the x— axis 



and the coupling interaction F is derived from the Lienard-Wiechert potential JT(J in the 
projectile frame 

Z T e~f 



<j>(x,y,z,t) 



(3) 



were Z^e is the charge of the target, h is the impact parameter, and 7 is the Lorentz factor, 



7 = - (v/cy. 

To be specific, we study the 208 Pb + 208 Pb collision at 640 A-MeV. We first solve the 
Schrodinger equation for the intrinsic motion, described by Hq. For this purpose it is con- 
venient to recast the intrinsic Hamiltonian into the following equivalent form 



Hn = huj 



^(7r 2 + p 2 )+/V 



(4) 



In the above, the commonly used variable transformations 

Pi 



Pi 



7Ti 



(5) 



V n " ' VdKJ 

have been made, where ti and pi stand for the components of the position and momentum 



operators respectively. The oscillator frequency is given by 



hu = hJ — 



and the dimensionless strength (3 is related to B as 

'A (huf D 2 



(6) 



B 



h 4 



(3. 



(7) 



In Fig. 1, we show the ratios E 1 ^q DR / (2 Egdr) an d E 1 ^q DR / (2 E GDR ) as a function of 
B, in the same range as chosen in Ref. 0. In this range, the anharmonicity can be treated 
using first order perturbation theory to great accuracy (~ 2%). The GDR and DGDR 
energies, to first order in f3, are given by 



E, 



(0) = hu (1 + 5/?) 



GDR 



E l D =g DR ([3)=2huj (1 + 7.5/5) 
E l ^ D M = 2hu J (1 + 6/3). 



(8) 
(9) 
(10) 



Fig. 1 is equivalent to that shown in Ref. and our results are essentially identical to 
theirs. 

The reduced transition matrix elements can also be easily calculated to first order in the 
parameter (3. We find 



a x 1/2 



(GDR \\E1 1| GS) = e i-^-j (1 - 2.5(3) , 
(DGDR, I = ||M|| GDR) = e \ J- (1 - 5(3), 

/ 9i \ 1 i 2 /To 

(DGDR, I = 2 \\E1\\ GDR) = e -1 J— (1 - 3.5/3) , 



where e is the absolute value of the electron charge and Si is given by the energy-weighted 
sum rule, 

S - 9 fr 2 NZ 

1 ~ 47r 2m A 

The energy-weighted sum rule for transitions from the ground state and from the GDR are 
satisfied to first order in the parameter (3, using the above energies and reduced matrix 
elements. 

In order to maintain Egdr{(3) at the experimental value, namely Egdr{P) = -E'gdr (13.4 
MeV, in the present case), the oscillator frequency must be renormalized as (3 is changed. 
The resulting renormalized frequency, from Eq.(|), is 

Note that in the 5-range of Fig. 1, the dimensionless parameter varies in the range —0.014 < 
(3 < 0.014 which yields 1.08 E C q]^ r > %uj{(3) > 0.93 E^ R . Whereas our oscillator frequency 
is a function of the anharmonicity parameter, in Ref. 0] it is kept constant at the harmonic 
value, Tiuj((3 = 0) = E^ R . This difference does not affect the ratio E l DGDR /{2 Eq DR ) shown 
in Fig. 1, since the oscillator frequency cancels out in this case (see eqs. (H) to (|T0D ). 
When the renormalized frequency is used in both the GDR and DGDR energies and matrix 
elements, the sum rules for transitions from the ground state and from the GDR are still 



satisfied. However, use of the renormalized frequency substantially changes the excitation 
probability of the DGDR, as will be shown below. 

The calculation of electromagnetic excitation probabilities and cross sections is performed 
with the code RELEX |9|], based on the Winther and Alder theory ||10|| . With this code, we 
perform a full coupled- channels calculation of the electromagnetic excitation of the GDR 
and DGDR. Similar results (about 10% larger) would be obtained when perturbation theory 
is used for the collision dynamics [fUJI . In Fig. 2, we show the enhancement of the DGDR 
excitation probability relative to its harmonic value as a function of B for the impact pa- 
rameter b = 30 fm. We find that for B ~ —100 MeV/fm 4 (which in this case corresponds to 
f3 ~ —0.7 xl0~ 2 ) the overall enhancement is 6%. For purposes of comparison, we have also 
performed calculations using a constant frequency (Hcu = 13.4 MeV in this case). We then 
obtain an enhancement of 35%, as shown by the dashed in line in Fig. 2, in agreement with 
Ref. (see their Fig. 1). 

In Fig. 3a, we show the enhancement in the impact-parameter integrated DGDR cross 
section (solid line) vs. B, for the same system. In the cross section calculations, impact pa- 
rameters up to 200 fm are taken into account and a lower cut-off at 15 fm is used to eliminate 
nuclear effects. The full line in Fig. 3a represents the result of the present work, in which 
an enhancement of only 4% is obtained for B = —100 MeV/fm 4 . The dashed line, obtained 
using a fixed value of the oscillator frequency, yields an enhancement of the DGDR cross 
section of 22% for the same value of B. The GDR cross section ratio cr GDR {B) / ctgdr{B = 0) 
obtained with fixed GDR energy, shown as a solid line in Fig. 3b, is close to one over the 
entire range of B values but is slightly less than one for large, negative anharmonicities, 
(about —0.5% at B = —100 MeV/fm 4 ). This small deviation is due to the increase in the 
population of the DGDR at these values of B and the correponding depopulation of the 
GDR. The GDR cross section ratio obtained with fixed oscillator frequency is shown as a 
dashed line in Fig. 3b. In this case, we find the GDR cross section to be enhanced by about 
10% at B = —100 MeV/fm 4 . The enhancement of 10% in the GDR cross section of Fig. 3b 
is clearly responsible for the large enhancement of 22% in the DGDR cross section of Fig. 
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3a at B = -100 MeV/fm 4 . 

The above conclusions do not change noticeably when the calculations are extended to 
other systems, such as 136 Xe + 208 Pb at 700A-MeV. The microscopic study of Ref. [12| 
established that the anharmonicity parameter scales as A~ l with the mass number. Thus, if 
B = —100 MeV/fm 4 represents a reasonable value for 208 Pb, then for 136 Xe a corresponding 
value would be B = —150 MeV/fm 4 . In Fig. 4, we display the results of calculations for this 
system as a function of the anharmonicity parameter B in Fig. 4. The solid line in the figure 
again shows the results of calculations in which the oscillator frequency is varied to maintain 
the GDR energy constant, while the dashed line represents the results of calculations in 
which the oscillator frequency is maintained fixed. Similar to the previous case, we find the 
enhancement of the DGDR cross section to be greatly reduced when the GDR resonance 
energy is maintained at a fixed value. As can be seen in Fig. 4, at B — —150 MeV/fm 4 , 
the DGDR cross section is enhanced by 62% when the oscillator frequency is maintained 
constant, but is enhanced by less than 10% when the GDR energy is maintained at its 
physical value. 

Before ending we comment briefly on the connection between the Bortignon-Dasso model 
used in this paper and microscopic models [OHT5J that aim to assess the importance of 



the anharmonic effects both on the spectrum and on the transition operator. Ref. fl4| 
finds, within the Lipkin model, small effects on the spectrum (which scale roughly as 1/A). 
Hamamoto finds, within nuclear field theory, that the nonlinear effects in the 1-phonon to 
2-phonon transition operator are also quite small and scale as 1/A ||16|| . As mentioned above, 
Ref. |fL2|| , through detailed microscopic calculations, finds that the anharmonic effects are 



indeed small and scale as 1/A. The values of the parameter B in both the Bortignon-Dasso 
and present calculations are taken to be small enough to be in line with the microscopic 
findings but also with the experimentally observed DGDR excitation energies (although the 
enhancement of the DGDR cross section could be increased thorough an artificially large B, 
there is no choice for this parameter that whould simultaneously explain the observed cross 
section enhancement and the only very small deviations of the DGDR excitation energy 
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from the harmonic limit). 

Another interesting point to mention is that the GDR has a width, which is considered 
neither by Bortignon and Dasso nor in the present calculation. The effect of the width of the 
GDR on the excitation of the DGDR has been recently studied within a harmonic picture 
[I7J. The overall effect of the width, at the energies considered here is, to produce a slight 



increase in the DGDR cross section, although not enough to explain all the available data. It 
would certainly be of interest to extend the present calculation within the anharmonic model 
by coupling the oscillator to other degrees of freedom (which would generate the damping 
width). 

In conclusion, we have investigated the effect of anharmonicities in the excitation of 
the DGDR in relativistic heavy ion collisions, with the same macroscopic model used by 
Bortignon and Dasso We point out that variations of the anharmonicity strength must be 
accompanied by a renormalization of the oscillator frequency, in order to maintain the GDR 
energy at a value consistent with the experimental one. We have found that this condition 
strongly reduces the enhancement in the DGDR excitation probabilities and corresponding 
cross sections, so that they remain much below the experimental results. 
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41.96.0886.00. D.T.P and L.F.C. acknowledge partial support from the Fundagao Uni- 
versitaria Jose Bonifacio, and M.S.H. and B.V.C. acknowledge support from the FAPESP. 
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Figure Captions 

• Figure 1: The ratio E l DGDR / (2 E GDR ) vs the anharmonicity parameter B, for 208 Pb. 
The solid line is for / = 2 and the dashed line for I = 0. The reduced mass for the 
oscillation of protons against neutrons is used for the mass parameter D. 

• Figure 2: The enhancement in the excitation of the DGDR in the collision of 208 P5-|_ 208 
Pb at 640 A-MeV for the impact parameter b = 30 fm. The solid line represents the 
results of the present calculation while the dashed line corresponds to a constant 
oscillator frequency. 

• Figure 3: Enhancement factor of the (a) DGDR and (b) GDR cross sections in the 
collision of 208 Pb + 208 Pb at 640 A-MeV. The dashed lines correspond to the results 
obtained with fixed oscillator frequency, while the full lines correspond to a fixed Eq DR . 

• Figure 4: Enhancement factor of the DGDR cross section in the collision of 136 Xe+ 208 
Pb at 700A-MeV. The dashed line corresponds to the results obtained with fixed 
oscillator frequency, while the full line corresponds to a fixed Eqdr- 
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FIGURE 3 a 
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FIGURE 3b 
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